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^ I Abstract 

CN , In this article, we assume that there exist scalar D* D* , D^Dl, B*B* and B*B* 

{^JQ' molecular states, and study their masses using the QCD sum rules. The numerical 

^ , results indicate that the masses are about (250 — 500) MeV above the corresponding 

<^ ; D* - D* , D* - D*, B* - B* and B* - B* thresholds, the y(4140)_is unhkely a scalar 

, D*Dl molecular state. The scalar D*D*, D*D*, B*B* and B*B* molecular states 

^vq ' maybe not exist, while the scalar D'* D' , D'lD' ^, B'*B' and B'*^B' ^ molecular states 

maybe exist. 
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L—i! 1 Introduction 



> 



Recently the CDF Collaboration observed a narrow structure (which is denoted as the 
' y(4140) now) near the J/il^(i) threshold with statistical significance in excess of 3.8 standard 

. deviations in exclusive — > J/tljcl)K~^ decays produced in pp collisions at ^/s = 1.96 TeV 

[1]. The mass and width of the structure are measured to be 4143.0 it 2.9 it 1.2 MeV 



o 



and II.7I5 Q ^ 3.7 MeV, respectively. The narrow structure y(4140) is very similar to the 



5.0 

■ charmonium-like state y(3930) near the J/i/juj threshold [21 [3]. The mass and width of 

g\ '. the y(3930) are 3914.611^ ± 2.0 MeV and 34^8^ ± 5MeV, respectively [3]. 

There have been several explanations for the nature of the narrow structure y(4140), 
such as a D*D* molecular state (HIHEIEIEIEKTO], an exotic (J^*^ = 1"+) hybrid 
^ . charmonium [5], a cess tetraquark state [11], or the effect of the J/ipfp threshold [12] . 

^ I The mass is a fundamental parameter in describing a hadron, in order to identify 

the y(4140) as a scalar molecular state, we must prove that its mass lies in the region 
(4.1 — 4.2) GeV. In Ref.|13j. we assume that there exists a scalar D*D* molecular state 
in the J/ip4' invariant mass distribution, and study its mass using the QCD sum rules. 
The numerical result indicates that the mass is about My = (4.43 it 0.16) GeV, which is 
inconsistent with the experimental data. The D*D* is probably a virtual state and not 
related to the meson y(4140) [1^. In this article, we extend our previous work to study 
the D*D*, D*D*, B*B* and B*B* molecular states in a systematic way considering the 
SU{3) symmetry and the heavy quark symmetry. 

In the QCD sum rules, the operator product expansion is used to expand the time- 
ordered currents into a series of quark and gluon condensates which parameterize the 
long distance properties of the QCD vacuum. Based on the quark-hadron duality, we can 
obtain copious information about the hadronic parameters at the phenomenological side 

MM- 
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The article is arranged as follows: we derive the QCD sum rules for the scalar D*D*, 
D*D*, B*B* and B*B* molecular states in section 2; in section 3, numerical results and 
discussions; section 4 is reserved for conclusion. 

2 QCD sum rules for the D^Dl, B*B* and B^Bl molec- 

ular states 

In the following, we write down the two-point correlation functions li ji^{p) in the QCD 
sum rules, 

nj/,(p) = iy"d^xe^P-^(0|T{jMx)J/,?t(o)}|o), (1) 

J(x) = Q{x)^^s{x)-s{x)rQ{x), (2) 
T]{x) = Q{xh^,u{x)dixh>'Q{x) , (3) 

where Q = c,b. We choose the scalar currents J{x) and r]{x) to interpolate the molecular 
states D*D*, D*D*, B*B* and B*B*, respectively. 

We can insert a complete set of intermediate hadronic states with the same quantum 
numbers as the current operators J(x) and 'r](x) into the correlation functions nj/^(p) to 
obtain the hadronic representation |14l [TS] . After isolating the ground state contributions 
from the pole terms of the scalar molecular states Y (we use the Y to denote the scalar 
molecular states D*D*, D*D*, B*B* and B*B*), we get the following result, 

n./.(p) = ' (4) 

where the pole residue (or coupling) Ay is defined by 

Xy = (o|JMo)|y(p)). (5) 

In the following, we briefly outline the operator product expansion for the correlation 
functions Ilj/j^{p) in perturbative QCD. The calculations are performed at the large space- 
like momentum region p"^ <C 0. We write down the "full" propagators Sij{x) and Cij{x) 
of a massive quark in the presence of the vacuum condensates firstly [15j . 

CI / \ ^^ij ^ijf^s ^ij /- \ , ^^ij I- \ / ^ijX ,_ ^ \ 

oijix) = — -TT-^ -(ss) -\ -ms(ss) X — (sOsCrGs) 

■ r- 2 

1152 ^ ^ 327r2x2 ''"^'^ ' ^ ' 

^ / X _ i f jiu^-ik-x I ^ij 9sGif + rnq) + (/f + mQ)a^^ 

+—(^^)^v^Q7i:2 — + (7) 



3^ 7T '-'^-^{k^ 
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where (sgguGs) = {sQsCTapG^^ s) and ( ) = {-^—^ ), then contract the quark fields 

in the correlation function Ilj{p) with Wick theorem, and obtain the result: 

nj(p) = iy"d^xe^P-"Tr[7^5i,(x)7aC,i(-x)]Tr[7'^C™„(x)7-5n,„(-x)] , (8) 

where the i, j, m and n are color indexes. Substitute the full s and Q quark propagators 
into the correlation function Ilj{p) and complete the integral in the coordinate space, then 
integrate over the variables in the momentum space, we can obtain the correlation function 
Ilj{p) at the level of the quark-gluon degrees of freedom. The correlation function n^(p) 
is calculated in the same way, we prefer neglect the technical details. 

In the QCD sum rules for the tetraquark states (irrespective of the diquark-antidiquark 
type and the molecule type) which have one or two heavy quarks, we always calculate the 
light quark parts of the correlation functions in the coordinate-space where the masses 
of the u, d, s quarks are taken as small quantities and treated perturbatively, and use 
the momentum-space expression for the heavy quark propagators |15j . then transform 
the resulting light-quark parts to the momentum-space with L'-dimensional Fourier trans- 
form [161 \T7\ Il8j . The main uncertainties in the QCD calculations originate from the 
high dimensional vacuum condensates, which are known poorly compared with the low 
dimensional vacuum condensates, for example, the quark condensate (qq) and the gluon 
condensate 

In this article, we carry out the operator product expansion to the vacuum condensates 
adding up to dimension-10 and take the assumption of vacuum saturation for the high 
dimensional vacuum condensates, they are always factorized to lower condensates with 
vacuum saturation in the QCD sum rules, and factorization works well in large limit. 
In the real world, Nc = 3, there are deviations from the factorable formula, we can 
introduce a factor k to parameterize the deviations, for example, 

(ss)"^ , {ss){sgsaGs) , {sgsCrGs)"^ Kiss)"^ , Hi{ss){sgs(jGs) , KisgsCrGs)"^ . (9) 

In Ref.[T^, we study the mass M^,^, with variation of the parameter k at the interval k = 
0—2. At the range = (2.6—3.0) GeV^, the value k = 1±1 leads to an uncertainty about 
50MeV, which is too small to smear the discrepancy between the theoretical prediction 
and the experimental data. If we assume the k has the typical uncertainty of the QCD sum 
rules, say about 30%, the correction is rather mild, we can neglect the uncertainty safely 
and take k = 1, i.e. the factorization works well. In the QCD sum rules for the masses of 
the p meson and the nucleon, the value of the k is always larger than 1 [19]. In calculations, 
we observe that larger k means slower convergence in the operator product expansion, 
requires larger threshold parameters, and results in larger ground state masses. We can 
draw the conclusion tentatively that the uncertainties (i.e. k > 1) in the QCD calculations 
enlarge the discrepancy between the theoretical prediction and the experimental data, our 
predictions based on the value k = 1 are reasonable. 

The contributions from the gluon condensates are suppressed by large denominators 
and would not play any significant roles for the light tetraquark states [20i |21] , the heavy 
tetraquark state [18] and the heavy molecular state |13j . In this article, we take into 
account them for completeness although their contributions are rather small. 

Once analytical results are obtained, then we can take the quark-hadron duality and 
perform Borel transform with respect to the variable = —p^, finally we obtain the 
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following two sum rules for the interpolating current J{x): 

Aye aP" = / dsp{s)e , 

J 4(mQ+ms)2 

where 



Pis) 



Pois) + P(ss) is) + /9(GG) is) + P{GG) is) { 



IT 



P{ss}ns) 

„B 



(10) 



(11) 



the lengthy expressions of the spectral densities /Oo(s), /0(ss)(s), P{GG)i^^'' P{GG)i^^ 
P{ss)'^is) are presented in the appendix. With a simple replacement, 



vfis , (ss) , (sgsCrGs) 



rriq , (qq) , {qgs<yGq) , 



(12) 



we can obtain the corresponding two sum rules for the scalar current r](x). 

Differentiating Eq.(lO) with respect to -p^, then eliminate the pole residues Ay, we 
can obtain two sum rules for the masses of the molecular states Y, 



so 



ds- 



4(mQ+m^)2 "'*d(-l/M2) 



:pis)e~ 



iSruQ+m.)^ dspis)e 



(13) 



the corresponding two sum rules for the scalar current r]{x) can be obtained analogously. 



3 Numerical results and discussions 

The input parameters are taken to be the standard values (qq) = —(0.24 it 0.01 GeV)^, 
(ss) = (0.8±0.2)(gg), (sgsaGs) = mg(ss), mg = (0.8 ± 0.2) GeV^, (^^) = (0.33 GeV^ 
niq ^0,ms = (0.14 ± 0.01) GeV, = (1.35 ±0.10) GeV and rub = (4r7±0.1)GeV at the 
energy scale /x = 1 GeV [12 ESI [22] . 

In the conventional QCD sum rules [14^ I15j . there are two criteria (pole dominance 
and convergence of the operator product expansion) for choosing the Borel parameter 
and threshold parameter sq. We impose the two criteria on the molecular states to choose 
the Borel parameter Af^ and threshold parameter sq- 

The contributions from the different terms in the operator product expansion are shown 
in Figs. 1-4, where (and thereafter) we use the (ss) to denote the quark condensates (qq), 
{ss) and the (sgsCrGs) to denote the mixed condensates (qgsCrGq), (sggCrGs). Prom the 
figures, we can see that the contributions from different terms in the operator product ex- 
pansion change quickly with variation of the Borel parameter at the values < 2.6 GeV^ 
and < 7.0 GeV^ in the hidden charm and hidden bottom channels respectively, such 
an unstable behavior can not lead to stable sum rules, our numerical results confirm this 
conjecture, see Figs. 6-7. 

The dominant contributions come from the perturbative term and the {ss) + {sgsCrGs) 
term; and the interpolating currents contain more s quarks have better convergent be- 
havior. The contribution from the terms involving the gluon condensate ( "°^'"^ ) are very 
small, the gluon condensate plays a minor important role. The vacuum condensates of the 
high dimension {ss)'^+{ss){sggaGs) serve as a criterion for choosing the Borel parameter 

and threshold parameter sq. 
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Figure 1: The contributions from the different terms with variation of the Borel parameter 
in the operator product expansion for the current c'y^ud^'^c. The A, B, C, D, E and 
F correspond to the contributions from the perturbative term, (ss) + {sgsCrGs) term, 
(2^) term, {^)+{^) [{ss) + {sgsaGs) + {ss)^] term, {ss)^+{ss){sgsaGs) term 
and (sgsCrGs)'^ term, respectively. The notations q, /5, 7, A, p and r correspond to the 
threshold parameters sq = 21GeV2, 22GeV2, 23GeV2, 24GeV2, 25GeV2 and 26GeV2, 
respectively. Here we take the central values of the input parameters. 
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Figure 2: The contributions from the different terms with variation of the Borel parameter 
in the operator product expansion for the current c'y^ss^'^c. The A, B, C, D, E and 
F correspond to the contributions from the perturbative term, (ss) + {sgsCrGs) term, 
(2^) term, {^)+{^) [{ss) + {sgsaGs) + {ss)^] term, {ss)^+{ss){sgsaGs) term 
and (sgsCrGs)'^ term, respectively. The notations q, /5, 7, A, p and r correspond to the 
threshold parameters sq = 21GeV2, 22GeV2, 23GeV2, 24GeV2, 25GeV2 and 26GeV2, 
respectively. Here we take the central values of the input parameters. 
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Figure 3: The contributions from the different terms with variation of the Borel parameter 
in the operator product expansion for the current Iry^ud^'^b. The A, B, C, D, E and 
F correspond to the contributions from the perturbative term, (ss) + {sgsCrGs) term, 
(2^) term, {^)+{^) [{ss) + {sgsaGs) + {ss)^] term, {ss)^+{ss){sgsaGs) term 
and (sgsCrGs)'^ term, respectively. The notations a, /3, 7, A, p and r correspond to 
the threshold parameters sq = 132 GeV^, 134GeV^ 136 GeV^, 138GeV^ 140 GeV^ and 
142 GeV^, respectively. Here we take the central values of the input parameters. 
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Figure 4: The contributions from the different terms with variation of the Borel parameter 
in the operator product expansion for the current Iry^ss^'^b. The A, B, C, D, E and 
F correspond to the contributions from the perturbative term, (ss) + {sgsCrGs) term, 
(2^) term, {^)+{^) [{ss) + {sgsaGs) + {ss)^] term, {ss)^+{ss){sgsaGs) term 
and (sgsCrGs)'^ term, respectively. The notations a, /3, 7, A, p and r correspond to 
the threshold parameters sq = 132 GeV^, 134GeV^ 136 GeV^, 138GeV^ 140 GeV^ and 
142 GeV^, respectively. Here we take the central values of the input parameters. 
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At the values M^-^ > 2.6 GeV^ and sq > 23GeV^, the contributions from the high 
dimensional condensates (ss)^ + {ss){sgsaGs) are less than 15% (4%) in the c^y^j^ud'^^c 
{c^^j_ss^^c) channel; the contributions from the vacuum condensate of the highest di- 
mension (sgscrGs)'^ are less than 3% in all the hidden charm channels, we expect the 
operator product expansion is convergent in the hidden charm channels. At the values 
^min — 7.0 GeV^ and sq > 136 GeV^, the contributions from the high dimensional con- 
densates (ss)'^ + (ss) (sgsCrGs) are less than 11% (5%) in the bj^ud'y'^b (b^f^ss'y^b) channel; 
the contributions from the vacuum condensate of the highest dimension {sQsCfGs)'^ are less 
than 7% in all the hidden bottom channels, we expect the operator product expansion is 
convergent in the hidden bottom channels. 

In this article, we take the uniform Borel parameter M^^„, i.e. M"^-^ > 2.6 GeV^ and 
^min — 7-0 GeV^ in the hidden charm and hidden bottom channels, respectively. 

In Fig. 5, we show the contributions from the pole terms with variation of the Borel 
parameter and the threshold parameter. The pole contributions are larger than 51% 
(55%) at the value M^„^ < 3.0 GeV^ and sq > 23GeV2 (24GeV2) in the c-f^udjf'c 
(cj^ssjf^c) channel, and larger than 52% (55%) at the value M^„^ < 8.0 GeV^ and sq ^ 
136 GeV^ (138 GeV^) in the h^^ud'y^b {h^,ss-/^b) channel. A gain we take the uniform 
Borel parameter M^^^., i.e. M^„^ < 3.0 GeV^ and M^^^ < 8.0 GeV^ in the hidden charm 
and hidden bottom channels, respectively. 

In this article, the threshold parameters are taken as sq = (24±l)GeV2, (25±l)GeV2, 
(138 ± 2) GeV^, and (140 ± 2) GeV^ in the c-f^udjf'c, cjf.ssjf'c, b-f/.ud-ff'b, and bj^ssjf'b 
channels, respectively; the Borel parameters are taken as = (2.6 - 3.0) GeV^ and 
(7.0 — 8.0) GeV^ in the hidden charm and hidden bottom channels, respectively. In those 
regions, the two criteria of the QCD sum rules are full satisfied [T5]. 

Taking into account all uncertainties of the input parameters, finally we obtain the 
values of the masses and pole residues of the scalar molecular states Y, which are shown 
in Figs. 6-7 and Tables 1-2. 

From Tables 1-2, we can see that the uncertainties of the masses My are rather small 
(about 4% in the hidden charm channels and 2% in the hidden bottom channels) while the 
uncertainties of the pole residues Ay are rather large (about (18 — 22)%). The uncertainties 
of the input parameters {{qq), {ss), {sggCrGs), (qgsCrGq), m^, rUc and mi,) vary in the range 
(2 — 25)%, the uncertainties of the pole residues Ay are reasonable. We obtain the squared 
masses My through a fraction, the uncertainties in the numerator and denominator which 
originate from a given input parameter (for example, {ss), (sggCrGs)) cancel out with each 
other, and result in small net uncertainty. 

At the energy scale /i = IGeV, ^ 0.19 [23], if the perturbative 0{as) corrections 
to the perturbative term are companied with large numerical factors, 1 + ^(s,mQ)^, the 
contributions may be large. For example, we can make a crude estimation by multiply- 
ing the perturbative term with a numerical factor, say 1 + ^(s,mc)^ = 2, in the hidden 
charm channel, the mass Mj-),jj, decreases slightly, while the pole residue Aj)*^* increases 
remarkably. From Figs. 1-4, we can see that the main contributions come from the pertur- 
bative terms, the large corrections in the numerator and denominator cancel out with each 
other. In fact, the ^(s,mQ) are complicated functions of the energy s and the mass mg, 
such a crude estimation maybe underestimate the 0{as) corrections, the uncertainties 
originate from the 0{as) corrections maybe larger. 

In this article, we also neglect the contributions from the perturbative corrections 
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Figure 5: The contributions from the pole terms with variation of the Borel parame- 
ter M^. The A, B, C and D denote the C7^nd7^c, c'y^ss^'^c, b^^ud'y^^b, and b^^ss^^b 
channels, respectively. In the hidden charm channels, the notations q, /?, 7, A, p and 
r correspond to the threshold parameters sq = 21GeV2, 22GeV^ 23GeV^ 24GeV2, 
25 GeV^ and 26 GeV^ respectively ; while in the hidden bottom channels they correspond 
to the threshold parameters so = 132 GeV^ 134 GeV^ 136GeV^ 138 GeV^ 140 GeV^ 
and 142 GeV'^ respectively. 
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0{a^). Those perturbative corrections can be taken into account in the leading logarithmic 
approximations through anomalous dimension factors. After the Borel transform, the 
effects of those corrections are to multiply each term on the operator product expansion 
side by the factor, 

where the Tj/^ is the anomalous dimension of the scalar interpolating current J/rj(x), 
the To„ is the anomalous dimension of the local operator On{0) in the operator product 
expansion, 

T{jMx)J/7?t(0)} = Cn{x)OniO), (15) 

here the Cn{x) is the corresponding Wilson coefficient. 

We carry out the operator product expansion at a special energy scale, say = 1 GeV, 
and can not smear the scale dependence by evolving the operator product expansion side 
to the energy scale M through Eq.(14) as the anomalous dimension of the scalar current 
J/r]{x) is unknown. Furthermore, the anomalous dimensions of the high dimensional local 
operators have not been calculated yet, and their values are poorly known. In this article. 



we set the factor 



1, such an approximation maybe result in some 



scale dependence and weaken the prediction ability; further studies are stilled needed. 

The central value of the present prediction M^^fj, = (4.48ib0.17) GeV is slightly larger 
than our previous calculation Mj^*£,* = (4.43ib0.16)GeV |13j . In the present work, we take 
a slightly larger threshold parameter sq = (25ib 1) GeV^ rather than sq = (24ib 1) GeV^ to 
take into account the SU (3) breading effects and enhance the contribution from the pole 
term. From Table 1, we can see that the central values of the possible scalar molecular 
states are about (250 - 500) MeV above the corresponding D* - D*, D* - D*, B* - B*, 
B* - B* thresholds respectively [24], the D*D*, D*D*, B*B*, B*B* are probably virtual 
states. In the constituent quark models, the energy gap between the ground state and the 
first radial excited state is about 500 MeV. The central values listed in Table 1 are below 
the corresponding thresholds of the first radial excited meson pairs. The scalar D*D*, 
D*D*, B*B*, B*B* molecular states maybe not exist, while the scalar D'*D'* , D'*D'*g, 
B'*B'* and B'*B'l molecular states maybe exist. 

In Refs.[71[T3]. the same current c{x)^fj,s{x)s{x)^^^c{x) is used to interpolate the narrow 
structure y(4140), however, the conclusions are quite different. The discrepancy mainly 
originates from the high dimensional vacuum condensates, the vacuum condensates of 
dimension-9,10 and the gluon involved vacuum condensates of dimension larger than 4 
are neglected in Ref.fT]. Those condensates are counted as 1, 0{^), 0{^), 0{^) 
respectively, and the corresponding contributions are greatly enhanced at small M^, and 
result in rather bad convergent behavior in the operator product expansion, we have to 
choose larger Borel parameter M^, one can consult the contributions from the (sggCrGs)'^ 
term in Figs.l-F.2-F,3-F,4-F for example. If we neglect the terms concerning those high 
dimensional vacuum condensates and choose the input parameters (especially the value of 
the nic) as Ref.[7j, the experimental data can be reproduced. As a byproduct, we can see 
that the scale dependence of the QCD sum rules only weakens the prediction ability mildly. 
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However, we insist on taking into account the high dimensional vacuum condensates, as the 
interpolating current consists of a light quark-antiquark pair and a heavy quark-antiquark 
pair, one of the highest dimensional vacuum condensates is (ss)"^ x ( "°^'^ ). 

The c-quark mass appearing in the perturbative terms (see e.g. Eq.(18)) is usually 
taken to be the pole mass in the QCD sum rules, while the choice of the rric in the leading- 
order coefficients of the higher-dimensional terms is arbitrary [25] . The MS mass rric 
relates with the pole mass m through the relation 



mc(m^) = rh 



vr \ 'K / 



-1 



(16) 



where K depends on the flavor number n/. In this article, we take the approximation 
rric ~ ^ without the corrections for consistency. The value listed in the Particle 
Data Group is mc^m^) = 1.271q5i GeV I23|) it is reasonable to take the value rUc = 
mc(l GeV^) = (1.35 it 0.10) GeV in our works. In Ref.[13], we also present the result with 
smaller value nic = 1-3 GeV, which can move down the central value about 0.06 GeV. The 
central value My = 4.37 GeV is stih larger than the D*D* threshold about 150 MeV. 

We can interpolate the scalar molecular states which consist of the scalar, pseudoscalar, 
vector, axial-vector and tensor meson pairs with the quark currents Qqq'Q, Qi'y^qq' i^^Q , 
Q^^qq'j^Q, Ot'^Ts^^'T/xTsQ and Qcr^'^qq'a^uQ, respectively. Those molecule type inter- 
polating currents relate with the diquark-antidiquark type interpolating currents through 
Fierz reordering in both the Dirac spinor space and the color space. 
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where A*^ = y f-^j the A" with a = 1, 2, • • • ,8 are the Gell-Mann matrixes. The A"^ with 

^4 = 2, 5, 7 are anti-symmetric and the A'^ with 5 = 0, 1, 3, 4, 6, 8 are symmetric. 

We usually take the diquarks as the basic constituents following Jaffe and Wilczek 
|26l |27] to construct the tetraquark states with the diquark and antidiquark pairs. The 
diquarks have five Dirac tensor structures, scalar C75, pseudoscalar C, vector 6*7^75, axial 
vector C7^ and tensor Ca^j^u, where C is the charge conjunction matrix. The structures 
C7^ and Ca^y are symmetric, the structures C75, C and 6*7^75 are antisymmetric. The 
attractive interactions of one-gluon exchange favor formation of the diquarks in color 
antitriplet 3c, flavor antitriplet 3/ and spin singlet Ig [28^ I29j. 

Naively, we expect the scalar tetraquark states with the structures C75A'^ — 75CA'^ and 
— have the smallest masses. In Refs.[30l I31j. we study the scalar and vector hid- 
den charm and hidden bottom tetraquark states which consist of C75A'^ — 75CA^ type and 
C7^A^ — CA"^ (and C7^75A'^ — 75CA'^) type diquark pairs respectively in a systematic way; 
and observe that the masses of the vector tetraquark states are about (0.6 — 0.7) GeV larger 
than the corresponding ones of the scalar tetraquark states. Furthermore, we observe that 
the scalar tetraquark states with the structure C75A^ — 75CA'^ have much smaller masses 
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Figure 6: The masses of the scalar molecular states with variation of the Borel parameter 
M^. The A, B, C and D denote the c'y^ud^'^c, c^y^ss'y'^c, h^^ud'^'^b, and h^^ss^'^b channels, 
respectively. 

than the corresponding ones with the structure — |30pl . From Eq.(17), we 

draw the conclusion tentatively that the Q^^qq'^^Q-, Qt'^Ts'Z^'t^jTsQ, Qcr'^'^qq'a^^Q type 
molecular states may have smaller masses than the corresponding Qqq'Q, Qi'^^qq'i^^Q 
type molecular states. The conclusion is not robust enough, detailed analysis with the 
QCD sum rules is still needed. 

The LHCb is a dedicated h and c-physics precision experiment at the LHC (large 
hadron collider). The LHC will be the world's most copious source of the h hadrons, 
and a complete spectrum of the h hadrons will be available through gluon fusion. In 
proton-proton collisions at y/s = 14TeV, the bb cross section is expected to be ~ 500/x6 
producing 10^^ bb pairs in a standard year of running at the LHCb operational luminosity 
of 2 X 1032cm-2sec-i [32]. The scalar D'*D'* D'lD'l, B'*B'* and B'lB'l molecular states 
predicted in the present work may be observed at the LHCb, if they exist indeed. We can 
search for the hidden charm molecular states in the DD, D*D*, DgDs, D^D*, J/ipp, J/ip4', 
J/tpuJ, r/cvr, r]c'r], ■ ■ ■ invariant mass distributions and search for the scalar hidden bottom 
molecular states in the BB, B*B*, BgBg, B*B*, Tp, T(j), Tlv, r/fevr, • • • invariant mass 
distributions. Those decays maybe take place through final-state re-scattering precesses 
with exchanges of the intermediate mesons a, vr, p, D, D* , ••• in the t channels. 

The QCD sum rules is just a QCD-inspired model, we calculate the ground state mass 

^The results with the structure CX^ — will be presented elsewhere. 
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Figure 7: The pole residues of the scalar molecLilar states with variation of the Borel 
parameter M^. The A, B, C and D denote the c'y^ud'y^c, C7^ss7'^c, Iry^ud'y'^b, and 
bj^ss^'^b channels, respectively. 



molecular states 


masses 


thresholds [23] 




4.38 ±0.18 


4.014 


c7qSS7"c 


4.48 ±0.17 


4.224 


bjaudj^b 


11.14 ±0.19 


10.650 


b')aSS"f°'b 


11.24 ±0.18 


10.831 



Table 1: The masses (in unit of GeV) of the scalar molecular states. 



molecular states 


pole residues 




5.1 ± 1.1 




6.2 ± 1.1 


b'jaud'y'^b 


2.7 ±0.6 


bjaSSj'^b 


3.2 ±0.7 



Table 2: The pole residues (in unit of 10 ^ GeV^ and 10 ^ GeV^ for the hidden charm 
and bottom channels respectively) of the scalar molecular states. 
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by imposing the two criteria (pole dominance and convergence of the operator product 
expansion) of the QCD sum rules. In fact, we can take smaller threshold parameter sq and 
larger Borcl parameter to reproduce the experimental value of the y(4140) as a scalar 
D*D* molecular state by releasing the pole dominance condition. Wc usually consult 
the experimental data in choosing the Borel parameter and the threshold parameter 
Sq. The present experimental knowledge about the phenomenological hadronic spectral 
densities of the multiqiiark states (irrespective of the molecule type and the diquark- 
antidiquark type) is rather vague. More experimental data are still needed. 

4 Conclusion 

In this article, we assume that there exist the scalar D*D*, D*Dl, B*B* and B*B* 
molecular states, and study their masses using the QCD sum rules. Our predictions depend 
heavily on the two criteria (pole dominance and convergence of the operator product 
expansion) of the QCD sum rules. The numerical results indicate that the masses are 
about (250 - 500) MeV above the corresponding D* - D*, D* - D*, B* - B* and B* - B* 
thresholds, the y(4140) is unlikely a scalar D*D* molecular state. The scalar D*D*, 
D*D*, B*B*, B*B* molecular states maybe not exist, while the scalar D'*D'* D'*D'l, 
B'*B'* and B'*gB'l molecular states maybe exist, and may be observed at the LHCb. 

Appendix 

The spectral densities at the level of the quark-gluon degrees of freedom: 
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